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We report combined first-principle and tight-binding (TB) calculations to simulate the effects of 
chemical edge modifications on structural and electronic properties. The C-C bond lengths and 
bond angles near the GNR edge have considerable changes when edge carbon atoms are bounded to 
different atoms. By introducing a phenomenological hopping parameter t\ for nearest-neighboring 
hopping to represent various chemical edge modifications, we investigated the electronic structural 
changes of nanoribbons with different widths based on the tight-binding scheme. Theoretical results 
show that addends can change the band structures of armchair GNRs and even result in observable 
metal-to-insulator transition. 



Carbon nanotubes (CNTs) with chiral vectors (n, m) 
are metallic when (2n + m) or equivalently (n — m) is 
a multiple of 3. Armchair CNTs with chiral vectors 
(n,n) are always metallic while the zigzag CNTs (n, 0) 
are metallic only when n is a multiple of 3. In general, 
the electronic properties can be modified by attaching 
atoms or molecules along CNTs' sidewalls using chem- 
ical methods The recent development of graphene, 
a sheet of an unrolled CNT, has attracted a great deal 
of research attention^ Graphene is a two-dimensional 
one-atom-thick carbon layer with carbon atoms arranged 
in a honeycomb lattice. Experiments by using the me- 
chanical method and the epitaxial growth method show it 
is possible to make GNRs with various widths Sev- 
eral experiments have also been performed recently to 
investigate the transport properties of graphene.— The 
results show that graphene is an interesting conductor 
in which electrons move like massless and relativistic 
particles. Several groups show STM/STS graphene im- 
ages with well-defined hydrogen-terminated edges in an 
ultrahigh- vacuum (UHV) environment^^ 

An interesting question to investigate is whether a 
metallic CNT is still metallic when it is unwrapped into 
a graphene nanoribbon (GNR) i n i 12 i 13 For simplicity, we 
chose the same definition for GNRs as for CNTs in this 
paper. An unwrapped (n, m) CNT, for instance, is called 
a (n, m) GNR. Previous theoretical studies based on 
the tight-binding (TB) approximation have examined the 
electronic states and energy dispersion relations of GNRs 
with the assumption that hydrogen atoms are attached 
to the carbon atoms on the GNR edge. These calcu- 
lations show that the metallic or insulating feature in 
GNRs are different from that in CNTkJi*l£d£ Zigzag 
GNRs (n,n), which correspond to unwrapped armchair 
CNTs with zigzag edges, are still metallic (if the spin de- 
gree of freedom is not considered). The electronic struc- 
ture of armchair (n, 0) GNRs depends strongly on the 
nanoribbon width i 14 i 15 i 16 Armchair GNRs are metallic 
when n = 3p + 1 (p is an integer) , and otherwise they 
are insulating. That is to say, GNRs can be made either 
as metallic or as semiconducting materials by controlling 



their width or chirality. This remarkable characteristic is 
very useful in making graphene-based nanoscale devices. 
In addition, GNRs are much easier to manipulate than 
CNTs due to their flat structure, and thus they can be 
tailored by using lithography and etching techniques. 

Because each edge carbon atom of GNRs is only 
bounded to two neighboring carbon atoms, a dangling 
carbon bond offers a remarkable opportunity to manip- 
ulate the electronic properties of GNRs. This can be 
done by attaching different atoms or functional molec- 
ular groups to the dangling carbon atom. Similar to 
the functionalization of CNT devices along edge dangling 
bonds j2a2*il the electronic properties of GNRs can also 
be alternated by chemical edge modifications. In this 
paper, we examine the geometric deformation of finite- 
width armchair GNRs caused by different edge addends. 
Our first-principle calculations show that the C-C bond 
length and bond angle near the edge undergo observable 
changes. To include the effect of the deformation on the 
graphene ribbon edge, we introduce a phenomenological 
hopping parameter t\ in our TB approximation calcula- 
tions. Our simulations show that the energy gap depends 
on ribbon width and hopping parameter t\. A nonzero 
energy gap exists for (3p +1,0) armchair GNRs, which 
means the metal-to-insulator transition can be achieved 
by edge modifications. 

The increments of C-C bonds and bonding angles at 
the nanoribbon edge have been reported based on the 
TB approximation calculations ! 11 ! 12 ! 13 . In this paper, we 
first estimate nanoribbon geometric deformation caused 
by various chemical addends. For a simple case (each 
armchair GNR edge carbon atom is saturated by one 
hydrogen atom), we evaluate geometric and electronic 
structure changes using the first-principle method. Our 
optimizations employ the Vienna ab initio simulation 
packag o 18 i 19 , which is implemented based on the local 
density approximation 20 of the density functional the- 
ory (DFT)*2i The electron-ion interaction is described 
by the ultra-soft pseudopotentials^ 2 - and the energy cut- 
off is set to be 286.6 eV. The atoms' positions are op- 
timized in order to reach the minimum energy with the 
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Hellmann-Feynman forces less than 0.02 eV/A. The re- 
sults show that the geometric relaxation localizes near 
the edge. Only the bond lengths and angles of edge car- 
bon atoms in armchair GNRs have considerable changes 
comparing with those of ideal graphene ribbons. For ex- 
ample, Z BAF=121.6° and ZBCD=118.3°. The inter- 
atomic distance between A and B sites (cIab) of (9,0) 
GNR is reduced from 1.42 A to 1.36 A, and d BC =1.40 A 
dcD = l-40 A and d£>#=1.42 A respectively, as shown in 
Fig.ffl 




Figure 1: Optimized structure of a (9,0) armchair GNR (all 
edge bonds are bound by hydrogen atoms). Here, the bond 
length is in A. Note that a zigzag (n, 0) GNR can be rolled to 
form an armchair ((n, 0) CNT, and an armchair (n, n) GNR 
can be rolled as a zigzag (n, n)) CNT. 

The geometric deformation of armchair GNRs, how- 
ever, depends on various kinds of chemical addends. For 
example, the edge carbon atoms connect to F atoms, the 
edge C-C bond lengths shorten to be 1.35 A (decreasing 
about 5%). In general, this kind of geometric deforma- 
tion results in the changes of hopping parameter between 
two neighboring carbon atoms on the GNR edge. The 
parameter change is defined as 



At 



< 2p / Jffi|2p' g > - < 2p z \H \2p z > 
<2p z \H \2p z > 



(1) 



Here, Hi and Hq are the Hamiltonians of the system 
with and without chemical edge modification, respec- 
tively. 2p' z and 2p z are the atomic orbitals of the cou- 
pled neighboring carbon atoms at the edge with the op- 
timized bond length and with the bond length of 1.42 A, 
respectively. Our calculations are conducted by using the 
SIESTA code in real space and single-zeta (SZ) basis* 2 ^ 
One can expect that the edge hopping parameter can in- 
crease (or decrease) as edge C-C distance is shortened (or 
expanded). The edge hopping parameter change for the 
hydrogen-saturated case is predicted to be 10.1 % based 
on the DFT, which is consistent with the numerical result 
based on the TB approximation* 2 ^ 

In previous TB approximation calculations ^ 12 ! 25 the 
dangling bonds are assumed to be saturated by hydrogen 
atoms and thus all transfer integrals between the nearest- 
neighbor sites are set to have the same values. This sim- 
ple choice of hopping parameter, however, does not con- 
sider geometric distortions at the nanoribbon edge. It 
is important to extend this existing scheme in order to 



understand the impacts of chemical edge modification on 
electronic properties of armchair GNRs. For simplicity, 
we adopt the TB approach to study these impacts. We 
choose the value of t\ either smaller or larger than the 
hopping parameter t = 2.66 eV of inner C-C bonds to 
simulate different chemical addends. Our theoretical re- 
sults show that only parameter t\ is sufficient to describe 
the electronic structure changes of armchair GNRs. 
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Figure 2: (color online) Top two valence bands and bottom 
two conduction bands of zigzag GNRs in t\ = t and t\ = 
1.2t cases. Here t = —2.66 eV is the hopping parameter of 
two-dimensional infinite graphene with the lattice constant 
CLc-c = 1.42. 



Next, we calculate the band structures of armchair 
GNRs, when edge GNR dangling bonds are bounded 
by atoms or molecular groups. We select ideal arm- 
chair GNRs with indices (n, 0) and calculate the band 
structures of three different ribbon widths (n = 8,9, and 
10) by setting t\ = t. The results are shown in Fig. 
2(a) to (c) for n = 8,9, and 10 with red solid lines, re- 
spectively. For clarity, we set the Fermi level to zero 
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(Ep = 0) and the wave number is normalized by the 
primitive translation vector for each GNR. In this study, 
we chose (8,0), (9,0) and (10,0) GNRs (their widths are 
around 2 nm) as examples to represent (3p + 2, 0), (3p, 0) 
and (3p+l, 0) GNRs (where p is an integer; and they cor- 
respond to three different width armchair GNRs labeled 
with (3p+l, 0), (3p, 0), and (3p+2, 0) in the literature 2 ^) . 
It is clear that the band structure of an armchair GNR 
depends on its width, and the energy gap is 0.58, 0.50, 
and 0.0 eV for (8,0), (9,0), and (10,0) armchair GNRs, 
respectively. Armchair GNRs (n, 0) are metallic when 
n = 3p + 1 as shown in Fig. 2(c). The energy gaps are 
nonzero for n = 3p and n = 3p + 2 GNR, in which GNRs 
become insulating. This result matches the previous re- 
ported results. 12 To simulate chemical edge modification 
effects, t\ is set to be 1.2t. The tt and ir* conduction 
bands are no longer degenerate at k = for a (10, 0) GNR 
and energy gap becomes 0.14 eV, which leads to metal-to- 
insulator transition. The DFT results suggest that both 
quantum confinement and edge effect cause the opening 
of energy gap. 26 Although (8, 0) and (9, 0) GNRs remain 
insulating when t\ = 1.2t, the energy gaps is changed to 
0.76 eV (increased) and 0.36 eV (decreased) by changing 
of the edge hopping parameter. Obviously, these theoret- 
ical results show that the electronic property of armchair 
GNRs is tunable via chemical edge modification. 

The band structure of zigzag (n, n) GNRs are similar 
to those of CNTs except for the existence of edge states 
that are caused by the gauge field at GNR edges i 12 i 13 
Due to the localized states at the edges, the uppermost 
valence band and the lowest conduction band are always 
degenerated at the Fermi level when fco < \k\ < tt. ko 
is slightly less than ^ for finite- width zigzag GNRs. If 
considering the spin degree of freedom in the DFT cal- 
culations, a hydrogen-saturated zigzag GNR is predicted 
to have a magnetic insulating ground state. 26 By select- 
ing different ti, we observe that the electronic structure 
near the Fermi level changes slightly, which indicates the 
electronic properties or carrier transport of zigzag GNRs 
are stable and insensitive to the edge modifications. 
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Figure 3: The energy gap of armchair GNRs depends on their 
widths, (a) when t\—t and (b) when t\ — 1.2t. 

Fig. 3 shows the relationship between armchair GNR 



widths (n) and their energy gaps when choosing t\ = t 
and t\ = 1.2t. For n = 3p, the energy gap of t\ = 1.2t 
is smaller than that of t\ = t. This value decreases when 
the width increases. For n = 3p + 1, the energy gap 
of t\ = t is zero, which becomes independent of ribbon 
width n. For t\ = 1.2t, C-C bonds on the edge are short- 
ened and the energy gap opens in (3p +1,0) armchair 
GNRs. For n = 3p + 2, the energy gaps of t\ = t are 
nearly the same as that of t\ = 1.2t. By comparing the 
electronic structure of t\ = t with those of t\ = 1.2t, we 
observe that the change of the edge hopping parameter 
(resulting from the added chemical groups) can signifi- 
cantly affect the electronic properties of armchair GNRs. 
Our calculations reproduce the results based on the local 
density approximation of the DFT. 26 It is worth men- 
tioning that the energy gaps of both t\ = t and t\ = \.2t 
approach zero, or the effect resulting from the addends 
becomes insignificant when n is very large. This obser- 
vation suggests that tuning band structure through edge 
modifications is effective only for finite-width GNR. 
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Figure 4: (color online) The effect of hopping parameter t± on 
energy gaps of armchair GNRs with different widths [ n = 8, 9 
and 10]. 

Next, we focus on evaluating the energy gap as a func- 
tion of hopping parameter t\ for three kinds of armchair 
GNRs ((8,0), (9,0),and (10,0)). Results are plotted in 
Fig. 4. It is interesting to note that around t\ = t (rep- 
resenting the slight deformation case), the energy gaps 
of (8,0), (9,0) and (10,0) armchair GNRs show differ- 
ent trends. In contrast to a (9, 0) GNR, the energy gap 
of a (8, 0) GNR increases as t\ increases around t\ = t. 
This trend is also appreciable in Fig.2. A (10,0) GNR 
has a zero energy gap at t\ = t i n i 12 i 13 Once we applied 
chemical edge modifications (t\ ^ £), there is an opened 
energy gap that always increases no matter whether t\ 
decreases {t\ < t) or increases {t\ > t). This result shows 
that a transition between metallic and insulating GNRs 
is achievable. Moreover, the energy gap can be controlled 
by selecting proper addends bounded to the edge carbon 
atoms of armchair GNRs. 

In Fig. 4, one t\ value exists that corresponds to the 
zero energy gap for both (8, 0) and (9, 0) armchair GNRs. 
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This phenomenon shows that, in principle, insulating 
armchair GNRs can be modified and become metallic 
GNRs. However, this kind of modification is difficult to 
implement since the hopping parameter t\ would have to 
be reduced or enlarged by 100%. 

In conclusion, first-principles calculations show that 
the chemical modification of armchair GNRs results in a 
considerable deformation of the bond lengths and bond 
angles near the edge. The introduction of hopping pa- 
rameter, ti 5 in the TB scheme accurately capture the 
effects caused by chemical edge modifications. Our 
theoretical results show that addends can change arm- 
chair GNRs' band structures and even lead to observable 
metal-to-insulator transitions. It should be pointed out 
that the chemical edge modification is effective only for 
finite-width GNRs. 

Note: After submitting this paper, we became aware 
of the similar independent work conducted by Son, Co- 
hen, and Louie et a/..— In their paper,— the scaling rules 
for GNR band gap as a function of their widths and their 



origins are studied based on the first-principles approach. 
A lattice model is then adopted within TB approxima- 
tion to explain the changes of energy gap in armchair 
GNRs by resolving the Hamiltonian perturbatively. Our 
numerical results presented in this paper are obtained by 
using the exact diagonalization technique. 

This work was partially supported by the National 
Natural Science Foundation of China with Grand num- 
bers 10574119, 10674121, 50121202, and 20533030, by 
National Key Basic Research Program under Grant No. 
2006CB9922000, by the USTC-HP HPC project, and by 
the SCCAS and Shanghai Supercomputer Center. Work 
at NTU is supported in part by COE-SUG grant (No. 
M58070001). Jie Chen would like to acknowledge the 
funding support from the Discovery program of Natural 
Sciences and Engineering Research Council of Canada 
under Grant No. 245680. The authors also would like 
to thank Thomas Mclntyre for the assistance with the 
finalization of this paper. 



Corresponding author. E-mail: liqun@ustc.edu.cn 
* Corresponding author. E-mail: jchen@ece.ualberta.ca 

1 C. W. Bauschlicher Jr., Chem. Phys. Lett. 322, 237 (2000). 

2 M. Burghard, Suf. Sci. Rep. 58, 1 (2005). 

3 X. Lu and Z. F. Chen, Chem. Rev. 105, 3643 (2005). 

4 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. 
Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, 
Science 306, 666 (2004). 

5 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
Nature 438, 197 (2005); Y. B. Zhang , Y. W. Tan , H. L. 
Stormer , Nature 438, 201 (2005); C. L. Kane , E. J. Mele 
, Phys. Rev. Lett. 95, 226801 (2005). 

6 Claire Berger, Zhimin Song, Xuebin Li, Xiaosong Wu, Nate 
Brown, Cecile Naud, Didier Mayou, Tianbo Li, Joanna 
Hass, Alexei N. Marchenkov, Edward H. Conrad, Phillip 
N. First and Walt A. de Heer, Science 312, 1191 (2006). 

7 C. Berger et al., J. Phys. Chem. B 108, 19912 (2004). 

8 Y. Niimi, T. Matsui, H. Kambara, K. Tagami, M. Tsukada, 
and H. Fukuyama, Appl. Surf. Sci. 241, 43 (2005). 

9 Y. Kobayashi, K. Fukui, T. Enoki, and K. Kusakabe, Phys. 
Rev. B 71, 193406 (2005). 

10 Y. Niimi, T. Matsui, H. Kambara, K. Tagami, M. Tsukada, 
and H. Fukuyama, Phys. Rev. B 73, 085421 (2006). 

11 H. Hosoya, H. Kumazaki, K. Chida, M. Ohuchi, Y. -D. 
Gao, Pure Appl. Chem. 62, 445 (1990). 

12 K. Nakada, M. Fujita, G. Dresselhaus, M. S. Dresselhaus, 
Phys. Rev. B 54, 17954 (1996). 

13 M. Fujita, K. Wakabayashi, K. Nakada, K. Kusakabe, J. 



Phys. Soc. Jpn 65, 1920 (1996). 

14 Y. -L. Lee, Y. -W. Lee, Phys. Rev. B 66, 245402 (2002). 

15 K. Wakabayashi, M. Fujita, H. Ajiki, M. Sigrist, Phys. Rev. 
B 59, 8271 (1999). 

16 Y. Miyamoto, K. Nakada, M. Fujita, Phys. Rev. B 59, 
9858 (1999). 

17 J. Taylor, M. Brandyge, and K. Stokbro, Phys. Rev. B 58, 
121101 (2003); M. Ernzerhof, M. Zhuang, and P. Roche- 
leau, J. Chem. Phys. 123, 134704 (2005). 

18 G. Kresse and J. Hafner, Phys. Rev. B 49, 14251 (1994). 

19 G. Kresse and J. Furthmiiler, Comput. Mater. Sci. 
6,15(1996); Phys. Rev. B 54, 11169 (1996). 

20 D. M. Ceperley and B. J. Alder,Phys. Rev. Lett.45, 
566(1980); J. P. Perdew and A. Zunger, Phys. Rev. B 23, 
5048(1981). 

21 P. Hohenberg and W. Kohn, Phys. Rev. 136, B 864 (1964); 
W. Kohn and L. J. Sham, ibid. 140, A1133 (1965). 

22 D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

23 J. M. Soler, E. Artacho, J. D. Gale, A. Garcia, J. Jun- 
quera, P. Ordejon, and D. Sanchez-Portal J. Phys.: Con- 
dens. Matter 14, 2745 (2002). 

24 D. Porezag, Th. Frauenheim, Th. Kohler, G. Seifert, and 
R. Kaschner, Phys. Rev. B 51, 12947 (1995). 

25 K. Sasaki, S. Murakami, R. Saito, Appl. Phys. Lett. 88, 
113110 (2006) 

26 Y. -W. Son, M. L. Cohen, and S. G. Louie, Phys. Rev. 
Lett. 97, 216803 (2006). 



